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Evolution of a system of interacting non-relativistic quantum flavor-mixed particles is consid- 
ered both theoretically and numerically. It was shown that collisions of mixed particles not only 
scatter them elastically, but can also change their mass eigenstates thus affecting particles' flavor 
composition and kinetic energy. The mass eigenstate conversions and elastic scattering are related 
but different processes, hence the conversion cross-section can be arbitrarily large even when the 
elastic scattering cross-section vanishes. The conversions are efficient when the mass eigenstates are 
well-separated in space but suppressed if their wave-packets overlap; the suppression is most severe 
for mass-degenerate eigenstates in flat space-time. The mass eigenstate conversions can lead to an 
interesting process, called 'quantum evaporation,' in which mixed particles, initially confined deep 
inside a gravitational potential well and scattering only off each other, can escape from it without 
extra energy supply leaving nothing behind in the well at t — > oo. Implications for the cosmic 
neutrino background and the two-component dark matter model are discussed and a prediction for 
the direct detection dark matter experiments is made. 



PACS numbers: 03.65.-w, 14.60.Pq, 14.80.-j, 95.35.+d 

I. INTRODUCTION 

A number of known and hypothetic particles are flavor- 
mixed, e.g., neutrinos, kaons, quarks, a neutralino, an 
axion (can be mixed with a photon), to name a few. 
How these particles behave in the non-relativistic limit 
has not been studied in detail. This paper addresses some 
important aspects of this profound physical problem. 

Mass (propagation) and flavor (interaction) eigenstates 
are the vectors obtained by diagonalizing the propagation 
and interaction parts of particle's Hamiltonian, respec- 
tively, and they can generally be not identical but related 
through a unitary transformation \fi) = \ m j)i 
where |/) and \m) denote the flavor and mass eigenstates, 
and U is a unitary matrix. Hence, a mixed particle pro- 
duced in a reaction has a specific flavor eigenstate, a, 
described by a wave-function being a superposition of 
several mass eigenstates pL. When a mixed particle is 
propagating, the mass eigenstates move with different ve- 
locities, which causes time-dependent interference known 
as flavor oscillations. 

An interesting and rather counter-intuitive property 
of non-relativistic flavor-mixed particles has been found, 
which is illustrated in the following example. Let us cre- 
ate a non-relativistic electron neutrino in a gravitational 
potential well. One should expect that if the neutrino is 
initially confined in the potential, it will remain confined 
forever (flavor oscillations do not change the picture). 
However, this is not so if the neutrino scatters elastically 
off other non-mixed particles from time to time. It has 
been shown that there is a non-vanishing probability to 
detect this electron neutrino outside the potential at a 
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later time, although no extra energy has been supplied 
to it 0. This effect, referred to as the "quantum evap- 
oration", is associated with mass eigenstate conversions 
— another process discussed in Ref. [2J. In our example 
here, a conversion of a heavier mass eigenstate yields a 
lighter one with a larger velocity. If this velocity exceeds 
the escape velocity, the light mass eigenstate is unbound 
and escapes to infinity. Note, however, that evapora- 
tion in such a thought experiment is not complete: only 
the heavy eigenstate can be converted into the escap- 
ing lighter eigenstate, whereas the initially created light 
eigenstate remains always bound if it was bound initially. 

In this paper, we consider a system of two non- 
relativistic flavor-mixed particles confined inside the 
gravitational potential well which can scatter off each 
other. We demonstrate that complete evaporation of 
both these particles is possible in this case. Moreover, 
we also show that the scattering and conversion cross- 
sections are fairly independent, so it is possible to have 
conversions even when the scattering cross-section van- 
ishes. We also show that the mass conversion cross- 
section in Minkowsky space (without gravity) is strongly 
suppressed in the mass-degenerate case. These results 
are important for better understanding of the proper- 
ties of mixed particle in general, as well as have interest- 
ing implications for the cosmic neutrino background and, 
possibly, dark matter physics and cosmology. 



II. INTERACTING MIXED PARTICLES 

In this paper we are interested in interactions of indi- 
vidual mass eigenstates. A mixed particle is created in 
a flavor state, but it consists of several mass eigenstates. 
Although all of these 'pieces' comprise a single particle, 
it is possible to visualize their kinematics as if they were 
normal particles having different masses. In the non- 
relativistic limit, these mass eigenstates propagate with 
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different velocities and separate from each other rapidly. 
(This is quite different form the relativistic case in which 
all eigenstates propagate nearly at the speed of light and 
it may take enormous time for them to separate, hence 
the plane wave approximation is commonly used.) There- 
fore, fairly soon, a mixed particle becomes a collection 
of spatially separated mass eigenstates which can inter- 
act with other particles independently. Apparently, an 
ensemble of non-relativistic flavor-mixed particles is, in 
most cases, an ensemble of individual mass eigenstates. 
Therefore, it is very natural to investigate the evolu- 
tion of such an ensemble in the mass basis rather than 
the flavor basis, which is usually used. The interaction 
matrix, however, is non-diagonal in the mass basis and 
off-diagonal terms represent transitions between different 
mass eigenstates. Should the mass eigenstates overlap 
to represent a particular flavor, these off-diagonal cou- 
plings 'balance' transitions of mass eigenstates into each 
other precisely to produce the scattered particle in a fla- 
vor eigenstate again. If, however, an individual mass 
eigenstate interacts, there is no such a 'balance', so new 
(absent) mass eigenstates are produced. Thus, one mass 
eigenstate can be converted into others. Such a process 
is of primary interest to us. 

Here we consider a simple model of stable two-flavor 
particles. We will interchangeably denote flavor eigen- 
states as \f a ) and \fp) or as just a and /?, whenever 
it's not confusing. Since masses of the mass eigenstates 
are different, we refer to them as heavy and light eigen- 
states, hence mi < m/,. Thus, similarly, the mass eigen- 
states are denoted as \rrih) and \mi) or as just h and I. 
A two-component flavor-mixed particle is described by a 
two-component wave-function, which representations in 
the flavor and mass bases are related via a 2 x 2 rotation 
matrix, U, where 9 is the mixing angle, i.e., 
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Fig. [T] illustrates such a particle. The bold red and 
blue curves represent heavy and light mass eigenstates 
assumed to have gaussian wave-packets, as in Eq. (26 1, 



and thin cyan and magenta curves are the corresponding 
flavor eigenstates; flavor oscillations occur where mass 
eigenstates overlap. 

Because each interaction involves two flavor-mixed par- 
ticles, the system is described by a two-particle wave- 
function, which has four components in the flavor and 
mass bases, namely 
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FIG. 1. Probability distributions of different eigenstates in 
space. The bold curves (red and blue) represent heavy and 
light mass eigenstates and thin curves (cyan and magenta) are 
flavor eigenstates. Flavor oscillations are seen in the overlap 
region. 



respectively, where the subscripts denote particle 1 and 
particle 2. Note that when the particles 1 and 2 are 
far apart (before or after an interaction), a two-particle 
wave-function is separable, being a direct product of one- 
particle ones: \f l fj) = |/<(xi,i)} ® |/j(x 2 ,f)), where % = 
a,/3, j = a, ft and \m p m q ) = \m p (s. u t)) ® |m g (x 2 ,t)), 
where p = h,l, q = h,l. The two-particle flavor and 
mass eigenstates are related as before, 
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in which c = cos 9 and s = sin 9. For simplicity, we 
will restrict further study to one-dimensional motion of 
particles. 

The evolution of the system at hand is described by the 
two-particle two-component Schrodinger equation. In 
the mass basis, it reads 



ihdt \mm(xx,%2,t)) — (H 
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Here the free particle Hamiltonian 
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satisfies energy conservation, where 
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and Am = rrih 



Jfgrav _ 



mi. Gravity enters via 
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where 
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m h (j)(xi) + m h (fi(x 2 ), 
m h (f)(xi) + mi4>(x 2 ) 1 
mi(f>(xi) + m h 4>(x 2 ), 
mi(j)(xi) + mi4>(x 2 ) 



and </>(a;) is an arbitrary gravitational potential. 

The interaction matrix is diagonal in the flavor basis, 
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where Vg a = V a p for indistinguishable particles. In the 
mass basis, we have 
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where the hermitian conjugate U 2 
real- valued unitary matrix and 
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Ul for the 



A=\ [3V aa + 2V a p + 3Vpp + 4(V4a - Vpp) cos 29 

+ {V aa -2V a p + Vpp) cos 49], 
B = § [W aa + 2V a0 + 3V p - 4(V aa - Vpp) cos 29 

+ {V aa -2V a p + Vpp) cos 49], 
C = § [V aa + 6V afj + Vpp - (V aa - 2V a p + Vpp) cos 49] 
D = \ [V aa - 2V a p + Vpp] sin 2 29, 



E=\ [V aa - Vpp + (V aa 
F = \ [V aa - Vpp - (V aa 



2^ + ^)008 20] sin 20, 
2V a p + Vpp) cos 29} sin26>, 



If the particles are distinguishable, one should make the 
substitution 2V a p — >• V a p + Vp a in the above equations. 
Since trace is invariant under a unitary transformation, 
Tr(V) = V aa + 2V a p + Vpp; also useful is . V? = 

vL + 2VZp + v$ p . 

The physics represented by the ^-matrix is easy to 
understand. There arc four different interaction combi- 
nations (input channels): hh —> . . . , hi — >• . . . , Ih — >• . . . , 
and 11 — > ... (although the states hi and Ih are identi- 
cal for indistinguishable particles by symmetry, we treat 
them separately for the sake of generality). In each 
of these interactions, there are four different outcomes 
(output channels): ••• — > hh, ■■■ — > hi, ■■■ — » Ih, 
and • • • — > 11. Thus, the V^-matrix 'sandwiched' be- 
tween initial and final states gives all 16 transition am- 
plitudes, Aki = i(mm\ V \mm) k , where k, I take the val- 
ues 1,2,3,4 = (hh),(hl),(lh),(ll). For example A 12 = 



(hh)(hl) 



(hl\V \hh) and corresponds to hh 



hi. 



These amplitudes are used to compute interaction cross- 
sections in the usual way [3]: for a given target particle 
t, the cross-section for a projectile particle of species s is 



'Siti-t-Sftf 
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where A is . ti){sftf) = (m Sf m tf \ V {sftf){Stt . } \m H m ti ) with 
Si and Sf being the initial and final states of the projec- 
tile particle, t% and t/ being those of the target particle, 
p Si and p Sf are the momenta of the projectile before and 
after the interaction, and are the elements of the inter- 
action potential matrix in the mass basis, Eq. (10). De- 



tailed calculation of matrix elements and cross-sections of 
mass eigenstate interactions is considered elsewhere [I]. 
Note that this form of a manifestly satisfies the principle 
of detailed balance [3], o^/p? = af^p-j, for the inte- 
gral cross-sections of forward and reverse processes, and 
that a oc 1/v. 

For completeness, we also present a useful parameter- 
ization of the cross-sections, which can easily be imple- 
mented in numerical models: 

Psr 

<J Sl u^s f t f = (7 — @{E S ft f ) B {StU)(sftf) , (12) 

where a is a common normalization parameter and 
Q(E Sftf ) is the Heaviside function which ensures that 
the process is kinematically allowed (i.e., negative final 
kinetic energy, E Sftf < 0, means the process cannot oc- 
cur), where E Sftf is equal to the initial kinetic energy of 
the particles in their center of mass frame plus AE, which 
is — Amc 2 for each I — > h conversion, or +Amc 2 for each 
h — > I conversion; thus, for example, for hi — > Ih, AE = 
and for in the process hh — > 11, AE = 2Amc 2 . The 
B-matrix with the elements B^i in the above equation 
are proportional to the squares of the properly normal- 
ized matrix elements of ^-matrix, Eq. (10 1, and simply 



enforces the types (or channels) of interaction that can 
occur. 

There are two types of processes: (i) elastic scatter- 
ings in which the system composition does not change 
(e.g., hh — ¥ hh, hi — ¥ Ih, etc.) and (ii) mass eigen- 
state conversions in which the composition changes (e.g., 
hh — > hi, 11 — > Ih, etc.). The diagonal elements or 
correspond to pure elastic scattering. Two off-diagonal 
elements, 23 and 32, describe 'mass exchange', but they 
contribute to scattering as well if the particles are indis- 
tinguishable. All other elements represent conversion of 
one or two mass eigenstates. The total energy and mo- 
mentum must be conserved in all processes. The energy- 
momentum conservation in elastic scattering is trivial, 
so we skip it. Conversions are different. Transitions in 
which a heavy eigenstate is converted into a light one 
go with the increase of kinetic energy and thus have no 
threshold. The opposite ones, where I is converted into 



h, have a threshold Amc 



mi)c and can only 



occur if kinematically allowed, i.e., if the initial kinetic 
energy of the interacting eigenstates is greater than the 
threshold. 
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Interestingly, there is a set of parameters, for which 
elastic cross-sections vanish identically but the conver- 
sion cross-section do not. Indeed, (i) the diagonal matrix 
elements, Eq. (10), namely A, B, C contribute to the 
total elastic scattering cross-section, cr scat ; (ii) the off- 
diagonal 'mass exchange' matrix elements V23 = V32 = D 
also contribute to er scat , if particles are indistinguishable; 
and (iii) the remaining elements E, F and V14 = V41 = D 
contribute to the total conversion cross-section, cr conv . It 
is easy to see that one can have cr sca i ~ simultane- 
ously with cr conv 7^ 0. First, scatterings like Ih — » Ih 
and hi — > ih vanish if C = D = 0, which requires that 
V a p — Vp a — 0, i.e., different flavors do not interact with 
each other, and also that Vg^ = — V aa . Second, scatter- 
ing channels hh — > hh and 11 — > 11 vanish if A = B = 0, 
which additionally requires maximal mixing, 9 = tt/4. 
Thus, the matrix V becomes 



V = V ac 



/o 1 1 0\ 

10 1 

10 1 

Vo 1 1 0/ 



(13) 



and V aa is the only independent matrix element. Thus, 
fscat = identically and <7 C onv oc |V aQ | 2 > 0, i.e., con- 
versions can occur even if the gas of mixed particles is, 
formally, collisionless. 



where we used that the velocities of h and I are also 
comparable, because ph = pi in the center of mass frame; 
hence v — p/m and v' = p'/m. Here we also intro- 
duced the "kick velocity", = c(2Am/m) 1 / 2 — this is 
the velocity a heavy eigenstate at rest gets upon conver- 
sion into a light eigenstate, provided recoil is vanishing. 
Thus, after the interaction, the mass eigenstates prop- 
agate along new geodesies, and if v[ is greater than the 
escape velocity of the potential, v[ > v csc , then both |m;)- 
eigenstates escape from the potential well. Alternatively, 
if elastic scattering occurs, hi — » hi or hi — > Ih, then 
the kinetic energy does not change and the eigenstates 
remain trapped. 

Therefore, upon any interaction involving the hi — > 11 
process the amplitude of the heavy eigenstate decreased 
irreversibly and both eigenstates can become unbound. 
Though the total probability remains unity, the probabil- 
ity to detect the particle (an electron neutrino, for exam- 
ple) inside the potential has decreased and the probabil- 
ity of its detection somewhere outside has become larger. 
Of course, the overall energy is conserved: the light eigen- 
state climbs up the potential and loses energy (e.g., a 
massless particle is redshifted). By repeating this cycle, 
one can further decrease the amplitude of the trapped 
eigenstate. Colloquially speaking, the particle (or parti- 
cles) "evaporates" from the potential well. 



III. KINEMATICS OF INTERACTIONS 

Let us consider an illustrative example of interaction 
of h and I belonging to two different particles. As we 
mentioned before, we consider one-dimensional motion, 
for simplicity. Let us consider hi — > 11 conversion and we 
assume here that the inverse process hi — » hh is kine- 
matically forbidden. Before the interaction, the mass 
eigenstates propagate along geodesies which are different 
(because the eigenstates have different velocities) and lo- 
calized in space (because the eigenstates are trapped). In 
the center of mass frame the momentum and energy con- 
servations axeph+pi = = p'i+p'i and (m 2 .c 4 +p 2 c 2 ) 1 / 2 + 
(mfc^+pfc 2 ) 1 / 2 = {mfc 4 +p[ 2 c 2 ) 1 / 2 + (rn 2 c 4 +p[ 2 c 2 ) 1 ' 2 , 
where 'prime' means 'after scattering', that is 

(ml^+pfc^ + irnfc'+pfcY/ 2 = 2(m 2 c 4 +pf c 2 ) 1 / 2 . 

and we remind here that the incoming particles are non- 
relativistic, Ph,Pi <C rrihC, mic. If rrih ^> mi, the outgoing 
mass eigenstates are relativistic 

^~£( ro 2+2m„m,-3m?)Va~!!|£. (15) 

Alternatively, if the masses are degenerate mi ~ ~ m 
and Am/m -C 1, then 

1/2 

Av = v' — v ~ [(Am/777,)c 2 + u 2 ] — v 

iiv«v k , 



IV. EVOLUITON OF A TWO-PARTICLE 
SYSTEM 

Here we show that two stable flavor-mixed particles, 
which are trapped in a gravitational potential and scat- 
ter off each other from time to time, gradually escape — 
or "evaporate" — from it. More precisely, the probability 
to detect the particles inside the potential decreases with 
time and the probability of their detection elsewhere in- 
creases. Such "evaporation" is a result of mass-eigenstate 
conversions in which a heavier eigenstate converts into a 
lighter one, thus adding kinetic energy to the scattered 
particles. We emphasize that the phenomena of mass 
eigenstate conversion and quantum evaporation are not 
related in any way to particle decays or other reactions, 
quantum tunneling and such. 

To illustrate the evaporation effect, we numerically 
solve the two-particle two-component Schrodinger equa- 
tion, Eq. (6j. To ease numerical computations, we 
chose a model potential with strong screening, <f>(x) = 
(f>oe~( x ' X9 > (1 — (x/xg) 2 ), where <^o < (meaning that 
the potential is attractive) determines its depth and x g 
sets its size (x g ~ 4 in computational units). The 
interaction potential is given by Eq. (10 1. Interac- 



tions of particles occur via a 5-function potential, i.e., 
Vaa, V a p, Vpp oc S(xi — X2), which is numerically repre- 
sented by V e-K* 1 -"*)/*"] 2 (1 + [{x l -x 2 )/x v ] 2 )- 1 , where 
Vb > and x v ~ 0.1; the actual shape of V{x\ — x<i) 
does not significantly affect the results so long ELS Xy IS 
small enough. The relative strengths are chosen to be 



5 




As a first case, we consider the interaction of a heavy 
and a light mass eigenstates belonging to two different 
particles, h\ + 1% — > . . . , which illustrates the effect of 
quantum evaporation. Fig. [2] shows the space-time dia- 
gram of the probability density of a heavy (orange) and a 
light (cyan) mass eigenstates given by Eqs. (18a), (18b I; 



FIG. 2. Space-time diagram of two interacting mass eigen- 
states in a gravitational potential well showing the effect 
of mass eigenstate conversions and particle evaporation, ob- 
tained by direct numerical solution of a two-particle two- 
component Schrodinger equation. The probability densities 
of light (cyan) and heavy (orange) mass eigenstates belong- 
ing to different particles are shown; yellow color originates 
from color blending of cyan and orange. The potential is lo- 
calized between x ~ —4 and +4 (in computation units). At 
t = the system consists of hi and h and both are trapped. 
During the evolution, each collision produces forward and re- 
flected wave-packets of all possible mass eigenstates; those 
corresponding to conversions escape to infinity in the form of 
light mass eigenstates. 



Vaa ■ V a p : Vpp =2:1: —2, the mixing angle is 
6 = ir/6 and the masses are chosen to be degenerate, 
Am/nih — 0.15. The initial wave-function components 
are taken to be gaussian wave-packets. 

Now we present exact numerical solutions of the 
Schrodinger equation for a pair of mixed particles. In 
order to simplify representation of the four-component 
two-dimensional time-dependent wave- function, Eq. (|3]), 
we compute probability densities of mass eigenstates for 
each particle (denoted by a subscript) as follows 

hh\{x,t) = j \hh(x, X2, t)\ 2 dx 2 , (17a) 

hh%{x, t) = / \hh(xi, x, t)\ 2 dx\ (17b) 



and similarly for other components. Since the particles 
are indistinguishable, we define the total probability den- 
sity of the heavy and the light mass eigenstates as 

h(x, t) = hhi + hh 2 + hh + lh 2 , (18a) 
l(x,t) = hl 2 + Ihi + lh + ll 2 . (18b) 



yellow color originates from color blending in regions 
where both mass eigenstates propagate along very similar 
paths. Initially, there is only the heavy mass eigenstate of 
particle 1 located at x = — 1 (in computational units) and 
the light mass eigenstate of particle 2 located at x = 1, 
and both are moving toward each other. Both eigenstates 
are initially gaussian wave-packets with momenta small 
enough to be trapped in the gravitational potential well. 
In each collision, forward and reflected wave-packets of 
all possible mass eigenstates are produced and light mass 
eigenstates participating in and/or resulting from conver- 
sions escape to infinity. 

To further elucidate the dynamics of the interactions, 
we show in Fig. ^ the wave-function components, 
Eqs. (17a), ( |17b[ ), namely, hhj (first panel), hlj (sec- 
ond panel) , Ihj (third panel) and 11 j (last panel), where 
j = 1,2. Here we use different color coding: orange 
represents particle 1 (i.e., j = 1) and blue represents 
particle 2. As one can see, at t = only the state hi 
(second panel) is non-vanishing; orange shows the wave- 
packet hl\ (i.e., the heavy eigenstate of particle 1 - the 
only heavy eigenstate initially present in the system) and 
blue is the wave-packet hl 2 (i.e., the light eigenstate of 
particle 2 - the only light eigenstate initially present in 
the system). The first interaction occurs at t ~ 0.2 and 
the second at t ~ 1.3. Note that several processes oc- 
cur at each interaction. First, no propagating hh wave- 
packets form, as is seen from the first panel, because such 
hi — > hh conversions are kinematically forbidden. Sec- 
ond, standard elastic collisions hi hi occur in which 
both forward- and back-scattered wave-packets are pro- 
duced, as is seen in the second panel. Third, elastic 
"exchange" hi — > lh also occurs, as is seen in the third 
panel (as we mentioned earlier, if the particles are indis- 
tinguishable, this process is equivalent to elastic scatter- 
ing). Here the wave-function Z/i-component, which was 
initially absent, appears at t ~ 0.2 as a vertex because 
both forward- and backward-scattered wave-packets ap- 
pear. After that, the wave-packets of both particles prop- 
agate but remain trapped in the potential, so they meet 
each other again at t ~ 1.1 (blue, particle 2) and t ~ 1.4 
(orange, particle 1). Note that although the wave-packet 
paths intersect, no interactions occur: the wave-packets 
belonging to the same particle do not self-interact but 
can only interfere. Finally, the fourth panel shows that 
light eigenstates are produced in conversions hi — > 11 seen 
as vertexes at t ~ 0.2 and t ~ 1.3. The velocities of these 
wave-packets exceed the escape velocity (controlled by 
the potential depth) so they leave the gravitational po- 
tential. 

Finally, Fig. [4] shows the expectation value of the num- 
ber of particles inside the gravitational potential. To sim- 
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FIG. 3. Space-time diagrams, as in Fig. but for the separate wave-function components: hh (left panel), hi (second panel), 
Ih (third panel) and 11 (right panel). Colors denote different particles: orange is particle 1 and blue is particle two. Initially, 
only hi component of the 2-particle wave-function is non-zero, meaning the presence of h mass eigenstate of particle 1 and I 
mass eigenstate of particle 2. Other mass eigenstates are produced upon interactions, which are seen as vertexes. See text for 
more details and discussion. 




FIG. 4. Normalized expectation values of the number of par- 
ticles confined inside the gravitational potential as a function 
of time, see Eqs. (19a|-(19cl. A heavy mass eigenstate is 



converted into a light mass eigenstate in each collision, at 
t ~ 0.2 — 0.3 and t ~ 1.3—1.4. Soon after each colli- 
sion/conversion, the light mass eigenstate escapes, thus de- 
creasing the total confined mass. 



plify comparison, we normalize them to the initial value 
as follows: 



h h (t) = 

hi(t) = 
h{t) = 



J h{x, t) dx 
J h(x, 0) dx ' 
J l(x, t) dx 
)l{x,Q)dx' 
hh + h u 



(19a) 

(19b) 
(19c) 



where Eqs. (|18aJ), (18b) were used. One sees that a light 
mass eigenstate is produced in each collision (t ~ 0.2—0.3 
and t ~ 1.3— 1.4) at the expense of the heavy eigenstate. 
Later, the light mass eigenstate escapes, thus decreasing 
the total mass inside. 

As a second case, we consider the full evolution of two 
flavor-mixed particles, each being a composition of both 
mass eigenstates. The essential difference of this case 
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FIG. 5. Space-time diagram of two interacting flavor-mixed 
particles in a gravitational potential, similar to Fig. [2] The 
probability densities of I (cyan) and h (yellow) mass eigen- 
states of both particles are shown. At t = the system con- 
sists of two particles of a particular flavor a, each being the 
superposition of hi, h and /12, h, and all the eigenstates 
are trapped inside a gravitational potential. At t > 0, each 
collision produces forward and reflected wave-packets of all 
possible mass eigenstates and those with v > v asc escape to 
infinity. 



from the previous one is that all mass eigenstates of both 
particles are present. The initial state of the system is 
two flavor-mixed particles produced as flavor eigenstates 
at x = — 1 and x = 1 for particle 1 and 2, respectively. 
These particles are a coherent mixture of mass eigen- 
states propagating with different velocities: the heavy 
eigenstates move toward each other and the light ones 
move initially away from each other. This initial setup 
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FIG. 6. Evolution of mass eigenstates shown which are overlaid in Fig. [5] Here, left panel (red) represents the h(x, t) eigenstate 
given by Eq. (18a| and right panel (blue) is the l(x,t), Eq. (18b I. 



allows us to separate the mass eigenstate interaction loca- 
tions thus simplifying the analysis of the dynamics. Fig. 
[5] is analogous to Fig. [2] and shows the conversion of 
heavy mass eigenstates into light ones and their escape 
from the gravitational potential. Cyan and yellow colors 
here denote I and h mass eigenstates. To elucidate the 
dynamics, we also separate the mass eigenstates into dif- 
ferent panels in Fig. [6j which is otherwise identical to 
Fig. Ill 



Figures [7J and [8] are similar to [3] and [4[ They show the 
evolution of the wave-function components and the num- 
ber of particles inside the potential. From Fig. [7] one can 
see that h mass eigenstates interact first (at t ~ 0.1, first 
panel) to produce elastically scattered trapped /i-states 
(panel one) and the outgoing h and Z-states via hh — > hi 
and hh — > Ih (panels two and three, respectively), both 
have large enough velocities to escape (h gets large v by 
recoil) . These escaping /i-states interact on their way out 
(at t ~ 0.4) with the scattered trapped /i-states to fur- 
ther produce escaping h- and Z-statcs via the processes 
of conversion hi — > 11, Ih — > 11 and "exchange" hi — > Ih, 
Ih — > hi; trapped /i-states are also produced at this time 
via inverse processes hi —> hh, Ih — > hh. Such processes 
repeat later as well, e.g., at t ~ 0.6. The amplitude of the 
direct hh —¥ 11 conversions is rather small for the chosen 
mixing angle and the values of ^-matrix, so they are not 
visible in this figure. However, they are seen in Fig. [8] 
as the decrease of the mass of the heavy eigenstate at 
t ~ 0.1, when only hh collision had occurred. The recoil 
velocity is larger in this process, hence the light eigen- 
state escape is fast. Overall, one can see from Fig. [8] 
that the particle evaporation is rapid and efficient in this 
case. 



V. ASYMPTOTIC STATE, t -> 00 

We demonstrated that evaporation of both light and 
heavy eigenstate can occur, which opens up a possibil- 
ity of complete evaporation of both particles, which were 
initially trapped. What conditions are needed for this to 
occur? Here we present some general estimates; a dedi- 
cated analysis may be needed for a specific system. Let 
the initial composition of the trapped particle popula- 
tion be rih,o an d n ;,o- For a single two-component parti- 
cle of flavor a, these are = c °s 2 & an d ni = sm2 ®i 
and for a particle of flavor /3, they are rih,o = sin 2 9 and 
nifl = cos 2 9, as follows from Eq. |lj. Note that in both 
cases rih,o + ni,o = 1> i- e -j there is exactly one particle 
in the system. If we consider a system of many parti- 
cles, rihfi and n^o must be multiplied by the number of 
particles. 

Let us also assume that the system is "optically thin" , 
i.e., probability of particle interaction during one bounce 
is very small, so if a conversion occurred, the escaping 
eigenstate experiences no further interactions and just 
leaves the system for good. We also assume that only 
forward conversions (h — > I) can occur; inverse processes 
(I — > h) are kinematically forbidden. We consider indis- 
tinguishable particles and also assume that Vk > v csc . 
These assumptions are very natural for non-relativistic 
mixed particles such as neutrinos (e.g., relic neutrinos 
from big bang) and some dark matter candidates because 
of their very small interaction cross-sections. 

The composition at t > is described by nu(t) and 
ni(t), which are governed by equations 

hh = -(&hh v ) n l - {°hiv) n h n u (20a) 
hi = -(<t m v) n h n u (20b) 

where we also assumed, for simplicity, that the particle 
density is uniform throughout the system. Here v is the 



FIG. 7. Space-time diagrams, as in Fig. |3| but for the two-particle case. Shown are wave-function components: hh (left 
panel), hi (second panel), Ih (third panel) and 11 (right panel); colors denote particles 1 (orange) and particle 2 (blue). Mass 
eigenstates conversions are seen as vertexes in these graphs (see text for detail). 




FIG. 8. Normalized expectation values of the number of par- 
ticles confined inside the gravitational potential as a function 
of time, similar to Fig. [4] Any time a heavy mass eigenstate 
is converted into a light mass eigenstate in a collision, soon 
after that the light mass eigenstate escapes to infinity. 



relative velocity of two interacting eigenstates which are 
comparable for heavy and light eigenstates if rrih — mi. 
Here also ath is the total cross-section of the processes 
hh — > hl,lh,ll and <7m is the total cross-section of the 
processes hi, Ih — > 11, hence 07^ oc 2E 2 + D 2 and 07,./ oc 
2F 2 , see Eqs. (10), (11). Whereas the general solution 



to these equations has no simple analytical solution, the 
asymptotic state can be found as follows. From Eqs. 



(20al, (20b) 



dn h _ a hh n h 
dm 



Phi ni 



+ 1. 



(21) 



This equation has a solution: 



nh(t) 

TlhQ 



1 - R 



1 



nt(t) 
ni, /n, h0 \ ( m(t) 



1-R 



ni,o 



(22a) 



where R = (J hh /a hl ^ 1, and 
n h {t) ni(t) 



n-h.o 



1 



"/.o , / ni(t) 



TlhQ 



(22b) 



if R = 1. We still do not know rih(t) and ni(i), but we 
note that h — > I conversions will occur until nh(t) 7^ 0. 
Therefore, asymptotically, when n^oo) — > 0, n;(oo) — > 
ni i00 - some constant value: 



nro 



(23a) 



which is valid for both < R < 1 and R > 1, and 



exp 



nh,o 
ni,o 



(23b) 



if R=l. 

We now conclude that when the initial composition 
satisfies the inequality 



nifi 
nh,o 



< 1 



&hh 

&hl ' 



(24) 



complete evaporation of mixed particles occurs, that is no 
particles will be left inside the gravitational well, n; ljOC = 
ni,ac = 0. Of course, the particles will be outside and 
traveling to infinity as light mass eigenstates only. This 
means that the flavor composition will be n a : np — 



sin 



cos 



VI. CONVERSIONS IN MINKOVSKY SPACE 

It is also important to investigate interactions of the 
particles in free space when gravity is negligible. This 
regime is relevant, for example, for the flavor-mixed dark 
matter in the early universe before structure formation 
starts, and for the relic cosmological neutrinos when they 
eventually become non-relativistic but still too hot to be 
confined by the gravitational attraction of the the large 
scale structure. 
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As before, mass eigenstates of a mixed particle move 
as if they are normal particles with certain (unequal) 
velocities and masses. The key difference between free 
and gravitationally confined particles is how their wave- 
packets spread with time. Depending on the shape of the 
potential, the wave-packet of a trapped particle, gener- 
ally, spreads slower than in free space or even contracts 
(e.g., near the turning points). In this case, the separa- 
tion of mass eigenstates occurs rapidly and can be nearly 
perfect as t — > oo, so one can treat these eigenstates in- 
dependently. In contrast, the wave-packets widths of free 
particles grow linearly with time and so does the sepa- 
ration between them. Therefore, the wave-packets of the 
two mass eigenstates can remain partially overlapped as 
t — v oo, and the effect may be very significant depend- 
ing on particle masses. Particle interactions in this case 
will involve both mass eigenstates leading to suppression 
of mass-conversion amplitudes. For example, then mass 
eigenstate wave-packets perfectly overlap, each particle 
is in a specific flavor eigenstate, and interactions do not 
change particle flavors (and hence mass eigenstate com- 
position) by definition of an eigenstate. 



Let us consider a non-relativistic mixed particle cre- 
ated at some moment of time t = at a position xq in a 
certain flavor eigenstate. It is a coherent superposition of 
mass eigenstates and each is described by a wave-packet, 
which we assume here to be gaussian: 



(27rAg)V4 



exp 



x ) 2 



4A 2 , 



(25) 



where A , rrij and Vj are the wave-packet width, mass 
and velocity and j — h,l. The first term describes a 
gaussian shape and the second term is simply the phase 
ikx = i(p/h)x — imvx/h. Note that Ao is the same 
for all mass eigenstates because the wave-packets must 
overlap completely at t = — the particle is created in 
a well-defined flavor eigenstate everywhere (i.e., at any 
x). Here we consider a one-dimensional case; general- 
ization to three dimensions is straightforward. At any 
time t > the wave-packet ipj(x,t) is given by the so- 
lution of the Schrodinger equation [5] for an initial state 
if>j(x) = i>j(x,Q), i.e., 



ipj(x,t) 



1/2 



2niht 
2ir ( A 



dx' exp 



rrij(x — x') 



l\2l 



2iH 



iht 



2toj ■ Aq 



-1/4 



exp 



(x — xq — Vjty 



rrijVj 



H — I rrijVjX J — £ -t 

4A§ + 2iht/m j h\ 3 2 



(26) 



This wave-function describes motion of j-th eigenstate 
with velocity Vj and the wave-packet spreading due to the 
momentum uncertainty, ApAo ~ h. In general, the ve- 
locities Vh , vi are different so the wave-packets of different 
mass eigenstates tend to separate in time: the gaussian 
centroids separate as 5x(t) ~ (vi — Vh)t oc t. On the 
other hand, the widths of the wave-packets also grow in 
time as A(t) ~ (h/rrijAo)t oc t as t — > oo. Since both 
grow linearly in time at late times, there will always be 
a non-zero overlap of the mass eigenstates. 

Interactions of mass eigenstates occur as follows. First, 
if the mass eigenstate wave-packets overlap completely, 
they both interact simultaneously as a flavor wave- 
function. This results in elastic scatterings only (flavor is 
conserved in interactions), because the interaction hamil- 
tonian, V, is diagonal is flavor basis, and no mass eigen- 
state conversions can occur. Second, in the opposite case 
of completely separated mass eigenstates, as in the case of 
trapping in a gravitational field discussed earlier, the in- 
teraction matrix is non-diagonal, so both elastic scatter- 
ing and conversions do occur. Finally, if the mass eigen- 
states partially overlap, there are non-zero chances for 
the particle to interact along both scenarios. In particu- 
lar, interactions as flavor eigenstates (i.e., non-separated 
mass eigenstates) is proportional to the overlap integral 



of the mass wave-packets. We calculate the overlap inte- 
gral now. 

A wave-packet given by Eq. ( 26 1 can be written as 
if)j(x,t) — Aj(x,t)e l 'l > *( 3 '' t \ where <j>j is a real- valued 
phase and A, is the real-valued amplitude which deter- 
mines the shape of the wave-packet. Since A| = ip*ipj 
we readily obtain: 



Aj = [27rA|(t)] 1/4 exp 



(x — Xq — Vjty 



where the wave-packet width is 



Af(i)=A 



2m j Aq 



4A?(t) 



if 



(27) 



(28) 



Note that xq and Ao are the same for both mass eigen- 
states because the particles are produced as flavor eigen- 
states, hence the mass eigenstate wave-packets com- 
pletely overlap at t = 0. If the particles form an ensemble 
in thermal equilibrium with some temperature T — the 
case that can be relevant to the early universe conditions 
— the expression for A(i) can readily be generalized [S] 
to yield 



A 2 Jt) = A 2 



kT 



<±m)Al 



(29) 
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The overlap integral of two mass eigenstates, h and I where 



I(t) 



A h (x,t)Ai(x,t) dx, 



(30) 



where Ah and Ai are given by Eq. (27). This integral is 
easily calculated analytically to yield: 



J(t) = 



A 2 



1/2 



exp 



t 2 fa-^) 2 

A 2 ) 



(AI 



(31) 



where A^ and A 2 are given by Eqs. (28) or (29). It's 
easy to check that 1(0) — 1, that is the wave-packets 
overlap completely at t — 0, and 1 > I(t) > at t > 0. 

To estimate the rate of mass eigenstate conversions, 
we look for the minimum overlap, i.e., for the asymptotic 
value of I(t) as t — > oo. We note that both mass eigen- 
states must have the same momenta p, so that = p/rrih 
and vi = p/mi. Note that if the momenta of the mass 
eigenstates are different, then their wave-packets, Eq. 
(26), carry extra ^-dependent phase even at t = 0, hence 



the particle exhibits flavor oscillations through space and 
this contradicts the physics of creation of the particle as 
a flavor eigenstate. 

First, we consider the case with strong mass- 
degeneracy: nth ~ mi Ri to, Am = mh — to; <C to. 
We obtain 



/(oo) 



1 



A TO 
TO 



Am" 



(32) 



where £ is a numerical factor of order unity. Indeed, if 
.Id: 

I p 2 A 2 



T = 0, Eq. (p8| holds, hence 

1 



p 



2h 2 



o(i), 



(33) 



where Ap is the momentum uncertainty and we used that 
ApAx ~ /i with Aa: ~ Aq and that Ap ~ p in collisions. 
In the opposite case when T is large enough for the first 
term in the brackets in Eq. (29) to be neglected, one has 



where — p 2 /2m ~ |fcT/m is the thermal energy of 
a particle. Overall, one can see that the value of /(oo) 
is fairly insensitive to the model assumptions and the 
estimate 



/(oo) ~ 1 - (Ato/to) 2 



(35) 



is robust. 

Second, if the masses are non-degenerate, m^ 3> mi, 
then 



\ m hJ 



1/2 / 3/2 s 



(36) 



/2 exp 



Ajp 2 
ft 2 



4A 2 fcTTOA 1/4 
ft 2 ) 



-y2e- 1 (l+4/3) 1/4 -G(l). (37) 

Thus, in this case the overlap is negligible, 

/(oo) ~ (to^to^) 1 / 2 « 1. (38) 

We have found that mass eigenstates can rapidly be- 
come well-separated in a gravitational field, where they 
propagate along significantly different geodesies, or in fiat 
space, where the local gravitational fields are extremely 
weak, provided there masses are very different. How- 
ever, if the mass eigenstates have degenerate masses and 
are propagating in flat space-time, their wave-packets 
spread much more rapidly than their centroids move 
apart. These mass eigenstates thus remain nearly per- 
fectly overlapped at all times, /(oo) ~ 1. Should it be 
identically unity, no conversions would occur. Due to 
the slight non-overlap, the conversion amplitude is small 
but nonzero, being a factor of (Ato/to) smaller than the 
conversion amplitude in the case of complete separation 
of the wave-packets. Thus the conversion cross-section in 
flat space, being proportional to the amplitude squared, 
is much smaller than that when mass eigenstates are well- 
separated, e.g., in the presence of sufficiently strong grav- 
itational field, thus 



if Ato -c to and a ' 



(ATO/TO) 4 CT conv , 

' fconv otherwise. 



(39) 



VII. IMPLICATIONS 

There are interesting cosmological implications of the 
obtained results. 

The first implication concerns with cosmological neu- 
trinos. Neutrinos from the cosmic neutrino background 
(CNB) have recently become non-relativistic; their ther- 
mal velocities are v t h — 81(1 + z)(eV/m v ) km s _1 [B], 
which is of the order of a few hundred to a thousand 
km/s, hence they can be trapped in dark matter halos 
of large galaxies and galaxy clusters. Scattering of neu- 
trinos off matter, though weak (but it can be greatly en- 
hanced by coherent effects [7]), will result in their mass 
eigenstate conversions and escape. The ultimate fate of 
the CNB will be that only the lightest mass eigenstates 
of neutrinos, \m) — (1,0,0), will exist in the universe as 
its age tji — y oo. From the Pontecorvo-Maki-Nakagawa- 
Sakata matrix, using presently measured mixing angles, 
we predict the asymptotic flavor composition, v e : : u T , 
to be 1 : (0.3 + 0.1 cos S) : (0.2 - 0.1 cos 5), where 5 is the 
yet-unknown CP-violating phase. Because of the very 
small neutrino cross-section and the very low average 
density in the universe, it will take much longer than the 
current age of the universe to achieve this asymptotic 
distribution. 
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Detectors on Earth, if they will ultimately be able 
to detect CNB neutrinos, should see the fractional de- 
viation from the uniform composition of order unity 
for upward vs. downward going relic neutrinos. In- 
deed, the non-relativistic neutrino-nucleon cross-section 
is CT ~ G 2 F El ~ 5 x 10- 56 (^/eV) 2 cm 2 with G F being 
the Fermi constant of weak interactions. Thus, for the 



heaviest species, assuming E 2 ~ Am?^ — 0.0027 eV , 
we have Co — 1-4 x 10 -58 cm 2 . The effect of coher- 
ent scattering increases the cross-section tremendously 
[7j: <t v ~ (JoZ 2 N 2 , where Z is the charge of atomic nu- 
clei, N ~ nV\ is the number of nuclei in the volume 
V\ ~ (47r/3)Aj B , n is the number density of nuclei and 
AdB ~ 0.5 cm is the neutrino de Broglie wavelength at 
z = (note, it is independent of m„ for CNB neutrinos). 
For Earth, Z ~ 25, n ~ 10 23 cm -3 , so the CNB neu- 
trino cross-section in Earth is ctcnb — 2 x 10 -10 cm 2 . 
The characteristic number density of the coherent scat- 
terers in Earth is n\ ~ 1/V\ and the typical distance 
neutrinos travel in Earth is its diameter, d ~ 10 9 cm, 
hence the 'optical depth' of Earth for the CNB neutri- 
nos is r ~ cfcNBn\d ~ 0.4, so the modifications to the 
relic neutrino composition and spectrum will be large. 
In particular, for the energy distribution of neutrinos one 
should expect spectral peaks at energies around yj Am 2 2 , 
yj Am 2 3 and yf Am 2 3 . 

The second implication is more speculative and deals 
with the recently proposed two-component dark matter 
[2 [8] . A number of dark matter candidates are flavor- 
mixed particles [§]. If two or more mass eigenstate are 
stable and haven't decay into the lightest one, also and 
these particles can self-interact, then the conversions dis- 
cussed in this paper will affect the composition, structure 
and dynamics of dark matter halos. The cold dark mat- 
ter (CDM) paradigm correctly describes the large scale 
structure of the universe but seems to fail at small scales. 
This is manifested by the departure of dark matter den- 
sity profiles in centers of halos (at scales less than tens of 
kiloparsecs in clusters and even smaller in galactic halos) 
from the 'cuspy' CDM profiles and the observed under- 
abundance of low- mass halos (with maximum circular ve- 
locities of less than a hundred kilometers per second) and 
the associated dwarf galaxies compared to the prediction 
of the CDM model. The weakly collisional two- (or multi- 
component dark matter (2cDM) model has a potential 
to resolve both problems. As in the self-interacting dark 
matter (SIDM) model JT], the central cusps are smeared 
out by collisions. However, unlike the SIDM model with 
a constant cross-section, the oc 1/v dependence of the 
flavor-mixed cross-section in the 2cDM model reduces the 
core size in galaxy clusters (which have an order of mag- 
nitude larger velocity dispersion than galactic halos) and 
brings them in agreement with observations. In addition, 
in the very centers of halos, where the 'optical depth' to 
collisions is much larger than unity, the dark matter will 
behave as a fluid and rapidly gravitationally collapse to 
form supermassive black holes |12j thus providing a pos- 
sible explaination for their existence at high redshifts. 



The observed paucity of small-mass halos can be due 
to the evaporation via h — > I conversions discussed 
in this paper. [Note, SIDM and similar models can- 
not address this problem at all.] This evaporation oc- 
curs if the characteristic velocity, Vk, the mass eigen- 
states get upon conversion is comparable or exceeds the 
escape velocity of the halo. The break in the max- 
imum circular velocity function of dark matter halos 
at v c ~ 50 — 100 km s _1 (corresponding to the halo 
mass of about 10 10 solar masses - this prediction of 
ours can be checked via cosmological observations) sug- 
gests the similar values for Vk being ~ v c , which im- 
plies the very high mass degeneracy Am/m ~ 10~ 8 . 
This is important for the following reason. The self- 
interaction cross-section should be rather large to be cos- 
mologically interesting: a s \jra > 0.1 cm 2 g _1 , but can- 
not exceed ~ 0(1) cm 2 g _1 to not contradict observa- 
tions. One would naively think that because of the large 
cross-section, all heavy eigenstates should be converted 
into the lightest one in the early universe soon after the 
dark matter freeze-out. This is not so, however, because 
the space-time is nearly flat in the early universe and the 
conversion cross-sections are strongly suppressed as com- 
pared to other interaction (e.g., scattering) cross-sections 
by a factor of (Am/m) 4 ~ 10~ 32 , thus making conver- 
sion processes irrelevant at that time. The possibility of 
strong mass degeneracy is interesting for the direct de- 
tection experiments as well. If the dark matter particle 
mass is ra ~ TeV, then the recoil energy detected in ex- 
periments can differ from the standard CDM prediction 
by ~ ±Am, that is by a few to 10 keV thus mimicing 
an inelastic or very light dark matter. It would be very 
interesting to look for such a signal in the ongoing and 
future experiments. 



VIII. CONCLUSIONS 

In this paper we studied the evolution of non- 
relativistic interacting flavor-mixed particles. We demon- 
strated that particle-particle interactions can lead to 
inter-conversions of their mass eigenstates, in addition 
to elastic scattering. These conversions are most effi- 
cient when the mass eigenstates are well-separated in 
space, but they are suppressed in flat space (no grav- 
ity) for mass-degenerate eigenstates, Eq. (39J). We stress 
that the conversions are not flavor oscillations: mass 
eigenstates remain intact during the oscillations. Also, 
conversions change momenta and kinetic energies of the 
eigenstates, Eqs. (15), (16), which results in the effect 
called "quantum evaporation" . Consider, for example, 
an ensemble of flavor-mixed particles trapped in a gravi- 
tational (or other non-flavor) potential well. Elastic colli- 
sions of the particles lead to mass eigenstate conversions, 
so the total number of particles trapped inside the well 
can decrease with time, Fig. [8j because of the escape 
of the conversion "secondaries". Whether the evapora- 
tion will be complete or else some particles will remain 
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trapped forever depends on the initial (flavor or mass) 
composition and the conversion cross-sections, Eq. (24 



We emphasize that unlike nuclear reactions, a partic 
kind does not change in conversions: neutrinos remain 
neutrinos, for example. What is changing is their flavor 
(and mass) composition, spatial localization, momenta 
and energies: some or all bound particles (i.e., with neg- 
ative kinetic plus potential energy) can become free (with 
positive energy) without extra energy supplied to the sys- 
tem. Obviously, evaporation is different from tunneling: 
particle's energy does not change in the latter. Finally, 
we discussed possible implications of the obtained results 



for (i) the cosmic neutrino background distortions in both 
the flavor composition and their energy and (ii) cosmol- 
ogy with two-component dark matter and possible reso- 
lutions of the core-cusp and substructure problems, and 
the better understanding of the early origin of supermas- 
sive black holes. A prediction for the direct detection 
dark matter experiments has also been made. 
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